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A new a-cluster wave function is proposed which is of the a-particle condensate type. Applications 
to 12 C and 16 show that states of low density close to the 3 resp. 4 a-particle threshold in both 
nuclei are possibly of this kind. It is conjectured that all self-conjugate An nuclei may show similar 
features. 
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Many properties of finite quantum systems such as nu- 
clei, quantum dots, atomic clusters or ultracould gases 
in a trap are fairly well described within a mean-field 
approximation neglecting the correlations between the 
quasiparticles. However, sometimes correlations become 
strong, giving rise to the formation of clusters, and they 
have to be taken into account. An intriguing problem 
is when bosonic clusters as bound states of fermions are 
produced, and the Bose character of the composite clus- 
ters competes with the fermionic properties of their con- 
stituents. As an example, we will discuss the relevance of 
a-like four-nucleon correlations in atomic nuclei. Special 
attention will be paid to such correlations which corre- 
spond to an a-type condensate in low-density symmetric 
nuclear matter, similar to the Bose-Einstein condensa- 
tion observed for finite numbers of bosonic atoms such 
as Rb or Na in traps. 

It is a well known fact that in light nuclei many states 
are of the cluster type jl] || . In the case of cluster states 
of stable nuclei where we have only very few excess nu- 
cleons in addition to the clusters, they are all located 
close to or above the threshold energy of breakup into 
constituent clusters. This fact which is known as the 
threshold rule |H) means that the inter-cluster binding is 
weak in cluster states. The threshold rule can be con- 
sidered as a necessary condition for the formation of the 
cluster structure, because if the inter-cluster binding is 
strong the clusters overlap strongly and the clusters will 
loose their identities. 

For example one of the fundamental questions of the 
cluster model is what kind of a-particle cluster states 
can be expected to exist around the threshold energy 
^na = n Ea of na breakup in self-conjugate An nuclei. 
One possible answer to this question, which is strongly 



under debate, is the existence of the cluster state of a 
linear na chain structure. The idea of the linear a chain 
state, originally due to Morinaga ||, is so fascinating 
that recently the formation of linear 6a chain states in 
24 Mg was studied extensively by experiments and also 
theoretically 0. The possibility of the linear 3a chain 
state in 12 C, which is the simplest linear a chain state 
was studied in detail by many authors solving the 3a 
problem microscopically However, these three-body 
studies all showed that the 3a-cluster states around the 
3a threshold energy Ef£ do not have a linear chain struc- 
ture. For example it was found that the calculated second 
+ state in 12 C, which corresponds to the observed sec- 
ond + state located at 0.39 MeV above the 3a threshold 
energy, has a structure where a-clusters interact predom- 
inantly in relative S-waves. Thus it was concluded that 
the cluster state near has not a linear chain struc- 
ture but rather an a-particle gas-like structure. 

On the other hand there have been recent theoretical 
investigations on the posibility of a-particle condensa- 
tion in low density nuclear matter In || Ropke et 
al. made a variational ansatz for the solution of the in- 
medium 4-body equation. In || Beyer et al. solved the 
Faddeev- Yakubovsky equations for an alpha-like cluster 
in nuclear matter. The outcome of these studies was that 
such a-condensation can occur only in the low-density re- 
gion below a fifth of the saturation value. At higher den- 
sities rather a state of ordinary p-n, n-n, or p-p Cooper 
pairing will prevail. In view of these results it may be a 
tempting idea that in finite self-conjugate An nuclei one 
could expect the existence of excited states of dilute den- 
sity composed of a weakly interacting gas of a-particles. 
Since the a-cluster is a Bose particle, such states could 
approximately be considered as an na cluster condensed 
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state and eventually excitations thereof. 

The purpose of this paper is to report on our study 
which not only confirms that indeed the second + state 
in 12 C could be considered as such a condensed state but 
that in addition also in 16 O such a state close to the 
threshold possibly exists. We will therefore then conjec- 
ture that the existence of such a-condensed states might 
be a general feature in N=Z nuclei. 

For the purpose of our study we write down a new type 
of a-cluster wave function describing an a-particle Bose 
condensed state: 

1$™} = (Ct)"|vac) (1) 
where the a-particle creation operator is given by 

Cl= J d 3 R e - R2 / R ° J d 3 n ■ ■ ■ d 3 r 4 

x <p Qs (r 1 -R)at iTi (ri)---^ 0s (r4-R)4 4 r 4 ( r 4) (2) 

with ^o s (r) = (l/(7r6 2 )) 3 / 4 e - r2 /(2fc 2 ) and ^(r) being 
the creation operator of a nucleon with spin-isospin ar at 
the spatial point r. The total na wave function therefore 
can be written as 

oc A{e-^^+- +x - > 0(ai) • • • <Ka„)}, (3) 

where B = (b 2 + 2i? 2 ) 1 / 2 and X, = (1/4) J2 n r ™ is the 
center-of-mass coordinate of the i-th a-cluster a{. The 
internal wave function of the a-cluster a^ is <p(oti) oc 
exp[-(l/86 2 ) Sl> n ( r ™ _ r m) 2 ]- The wave function of 
Eq. (3) is totally antisymmerized by the operator A. It is 
to be noted that the wave function of Eqs.(l,3) expresses 
the state where na-clusters occupy the same Os harmonic 
oscillator orbit exp[— -Jj-X 2 ] with B an indepedent vari- 
ational width parameter. For example if B is of the size 
of the whole nucleus whereas b remains more or less at 
the free a-particle value (a situation encountered below), 
then the wave function (3) describes an na cluster con- 
densed state in the macroscopic limit n — > oo. For finite 
systems we know from the pairing case that such a wave 
function still can more or less reflect Bose condensation 
properties. Of course the total center-of-mass motion can 
and must be separated out of the wave function of Eq.(l) 
for finite systems. In the limiting case of B = b (i.e. i?o 
= 0), Eq.(3) describes a Slater determinant of harmonic 
oscillator wave functions. It may also be worth mention- 
ing that for B ^ the wave function (1,3) is different 
from the cluster state proposed by Brink Q . 

The state |$ nQ ) has spin-parity + . In the limit of 
Rq = 0, the normalized wave function \$na) — \&na) 
/ v (®na\®na) is identical to a harmonic oscillator shell 
model wave function with the oscillator parameter b. 
For n — 3, it is identical to the p-shell wave function 
(0s) 4 (0p) 8 , [444] 0+), and for n = 4, it is identical to 



the double closed shell wave function, |(0s) 4 (0p) 12 , + ). 
This is easily proved by noticing that these limit wave 
functions for n =3 and 4 have maximum spatial sym- 
metry [444] and [4444], respectively. Only 8 Be has an 
a-particle structure in its ground state. Heavier na nu- 
clei collapse to the dense state in their ground state but 
the individual a's may reappear when these nuclei are 
dilated i.e. excited. 

We calculated the 

energy surfaces in the two parameter space, Rq and b, 
E na (R ,b) = {^ a (Ro,b)\H\^ a (R ,b)}, for n = 3 and 
4. The Hamiltonian H consists of the kinetic energy, the 
Coulomb energy, and the effective nuclear force named 
Fl which was proposed by one of the authors and con- 
tains a finite range three-nucleon force in addition to the 
finite range two-nucleon force jl0| . This force reproduces 
reasonably well the binding energy and radius of the a- 
particle, the a-a phase shifts of various partial waves, 
and the binding energy and density of nuclear matter. 
As we will see below this force also gives good results for 
binding energies and radii of 12 C and 16 O. 

In Fig.l we give the contour maps of the energy sur- 
faces E na (Ro,b) for 12 C and 16 0. The qualitative fea- 
tures of both surfaces are similar. They show a valley 
running from the outer region with large Rq > 1 1 fm and 
b k, b a — 1.44 fm to the inner region with small Rq and 
b > b a , where b a is the oscillator parameter of the free 
a-particle. The valleys have a saddle point at Rq « 10 
fm for n = 3 and at Rq ~ 10.6 fm for n — 4. Beyond the 
saddle point, E na (R ,b a ) » — nE a , where E a = - 
27.5 MeV is the theoretical binding energy of the free a 
particle by the present Fl force in Hartree-Fock approxi- 
mation. Therefore we have 3E a = -82.5 MeV and 4E a = 
-110 MeV. In the asymptotic region the average inter-a 
distance is large and the kinetic energy of the center-of- 
mass motion of an a- cluster (3^ /4mS 2 ) is very small 
which leads to E na (Ro,b a ) being more or less equal to 
E^. The height of the saddle point measured from the 
theoretical threshold energy is about 1.4 MeV for n = 3 
and 2.2 MeV for n = 4. The appearance of the saddle 
point is due to the increase of the Coulomb energy and 
kinetic energy towards the inward direction which is not 
yet compensated by the gain in potential energy around 
the saddle point region. This saddle point will help to 
stabilize the possible a condensed state around E™£ . The 
minimum of the energy surface is located at Rq ~ 2 fm 
for n = 3 and at Rq « 1 fm for n = 4. Since Rq = 
means the shell model limit, we thus see that the wave 
function even at the energy minimum point deviates from 
the shell model limit and shows rather strong a-particle 
correlations. The gain in energy from the shell model 
limit is 10.3 MeV for 12 C and 4.7 MeV for 16 0. Before 
comparing numbers with experiments we have to make a 
quantum mechanical calculation. This will be achieved 
via a standard Hill- Wheeler ansatz taking Rq and b as 
the Hill- Wheeler coordinates. However, in order to re- 
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duce the complexity of the calculation and because the 
valleys run essentially parallel to the Rq axis at b = b a 
we take b = b a = constant and only discretise the Rq 
variable. We therefore have 



X>((ik)i,M*™((*o)i>&«)>- 



(4) 



The normalization of fk({Ro,b)j) is so that the fc-th 
eigen- function |^,i Qi fc) is normalized. The adopted mesh 
size of Rq values is typically 0.5 fm. In order to see the 
character of the obtained wave function \^ n a.k), we in- 
troduce the overlap amplitude 



A na , k (R ,b) = (<f>% a (Ro,b)\y na , k ) 



(5) 



From this overlap amplitude we can estimate the relevant 
values of the variational parameter Rq in the different 
states k as will be discussed below. 

After outlining the results for the new kind of wave 
function for 12 C and 16 O, we will discuss whether the ob- 
tained condensed states correspond to the states found in 
these nuclei. We first consider 12 C, i.e. n = 3, see Table I. 
The calculated lowest two eigenenergies are situated at - 
85.9 and -82.0 MeV. The lowest energy state corresponds 
to the ground state of 12 C and is only slightly lower than 
the minimum point of the energy surface located at -85.5 
MeV. However, the calculated ground state is still above 
the observed binding energy of 12 C which is at -92.16 
MeV. An increase of mesh points will certainly lower the 
energy but, as has been discussed by many people, in 
order to reproduce the observed 12 C binding energy sat- 
isfactorily we have to extend our functional space so as 
to include the spatial symmetry broken wave functions 
which allow to incorporate the effect of the spin-orbit 
force adequately. The second eigenvalue lies 0.36 MeV 
above our theoretical 3a threshold energy, = -82.5 
MeV, and we believe that it corresponds to the observed 
second + state of 12 C which lies 0.5 MeV above Ef£ . 
As seen in Table I the rms radius of the obtained wave 
function ^30,2) is 4.29 fm which is much larger than the 
one of the ground state which is 2.97 fm, slightly greater 
than the experimental value 2.45 fm but in agreement 
with the missing binding of 6.75 MeV. We thus see that 
the second + state corresponds to a very dilute system 
of average density which is only about a fifth of the ex- 
perimental ground state density. 

To characterise the wave function by a typical value of 
the width parameter Rq, we consider the overlap ampli- 
tude Asa^iRo, b a ) given in Eq.(6) as a function of Rq at 
fixed b a . Whereas the ground state (k = 1) wave func- 
tion is almost exhausted by one ^^ (Rq, b a ) with Rq » 2 
fm which is quite close to the wave function of the mini- 
mum energy point of the energy surface, the second 0^ =2 
state has the largest overlap amplitude (about 0.87) with 
Rq w 4.5 fm. This rather large value implies that the dis- 
tribution of the center-of-mass momenta is rather narrow, 



in a certain approximation to an a condensate in infinite 
nuclear matter where all a-clusters populate the same 
state P = of the center-of-mass momentum. The fact 
that the calculated 0^ =2 state is of dilute density is in 
agreement with nuclear matter calculations |Bp| where 
it was shown that a condensate of a-like particles (quar- 
tering) is possible only in matter with density p < 0.03 
fm~ 3 . The average distance of the a-clusters in the dilute 
0fc=2 state is in agreement with this value for low den- 
sity nuclear matter, where the overlap of the a-clusters 
is small so that the Pauli blocking effects are weak. 

Let us now discuss the case of 16 0, i.e. n = 4. The 
energies of the lowest observed + states are shown in Ta- 
ble II, together with the corresponding widths. The first 
excited 0^ state at 6.06 MeV is very well known to have 
a-clustering character ]l],0] and is well described by the 
12 C + a microscopic cluster model as having the struc- 
ture where the a-cluster moves in a S state around the 



12 C-cluster in its ground state 
tcr states have been proposed M 



though also other clus- 
Similarly, the third ex- 
cited 0^ state at 12.05 MeV can be described by the same 
model where the a-cluster moves in a D state around the 
12 C-cluster in its first 2 + excited state pjj. We will ex- 
clude these well understood states from our further dis- 
cussion. The excited states 0^" at 11.26 MeV and 0^ state 
at 14.0 MeV observed in 12 C + a elastic scattering IllJ] 
cannot be described by such a model. Furthermore, they 
have very large decay widths, not typical for the other 
states. These states may be described by our new wave 
function as condensed states. 

As seen from Table I, the experimental value of the 
ground state (Of) of 16 at -127.62 MeV is well repro- 
duced by the calculated energy value for the ground state 
(0fe=i) & t -124.8 MeV. The calculated energy is above the 
minimum energy of the energy surface. It is because the 
b value of the minimum energy point is fairly larger than 
b a and the minimum energy point is not covered by the 
adopted mesh points. In order to have a better wave func- 
tion for the ground state we, of course, need to include 
mesh points around the minimum energy in our genera- 
tor coordinate calculation. When we adopt b = 1.57 fm 
which is the b value of the energy minimum point, the 
generator coordinate calculation gives -128.0 MeV as the 
lowest eigen energy. The rms radius of the calculated 
(0fe=i) state is 2.59 fm and is slightly smaller than the 
observed (0^~) rms radius, 2.73 fm, of 16 0. The second 
(0^_ 2 ) state of our calculation is bound by 6 MeV below 
the theoretical 4a threshold energy. The rms radius of 
this state is 3.12 fm and this state has the largest overlap 
amplitude (about 0.86) with $f a (R ,b a ) with Rq ps 2.5 
fm. We conjecture that this state corresponds to the ob- 
served O3" state situated at 3.18 MeV below the observed 
4a threshold energy. Indeed one may argue that there 
will be some mixing between the second (0i_ 2 ) state and 
the 12 C+a state, bringing theoretical and experimental 
energies closer together. The third (0jjL 3 ) state of our 
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calculation is bound by 0.7 MeV below the theoretical 
4a threshold energy. We think that it may correspond to 
the measured 0^" state situated at 0.44 MeV below the 
observed 4a threshold energy This state has a very large 
rms radius of 3.94 fm and has the largest overlap ampli- 
tude (about 0.80) with <f>% a (R , b a ) with R w 4.1 fm. In 
analogy to the case of 12 C, these values indicate that this 
state of dilute density should be considered as just the 
4a-cluster condensed state that we expected. One should 
point out that the ease with which we get the + -states 
around E££ is a strong indication that our wave func- 
tion (1,3) grasps the essential physics because otherwise 
the threshold states are very difficult to obtain. We also 
would like to mention that the present formalism yields 
very good results for the groundstate of 8 Be as well. 

In conclusion, our present study thus predicts in 12 C 
and 16 O the existence of near-na-threshold states which 
are the finite system analogues to a-cluster condensation 
in infinite matter. They are characterized by low density 
states so that the a-clusters are not strongly overlapping 
and by an n-fold occupation of their identical S-wave 
center-of-mass wave function. Therefore these states are 
quite similar in structure to the Bose-Einstein condensed 
states of bosonic atoms in magnetic traps where all atoms 
populate the same lowest S-wave quantum orbital. Be- 
cause of the short life-time of the a-condensed states, the 
predicted large values for the rms radii may be verified 
by indirect methods. The measurements of the spectra 
of the emitted a particles should allow to determine the 
Coulomb barrier which is expected to be small for the 
low density states. Of particular interest would be a-a 
coincidence measurement of decaying condensed states. 

We conjecture that such condensed a-cluster states 
near the na threshold may also occur in other heav- 
ier An self-conjugate nuclei. For example condensed 6a 
states of 24 Mg could be deformed and the measurement 
of a reduced moment of inertia over the rigid body value 
would be a strong indication for a-particle superfluidity 
The wave function we have proposed in this work is very 
flexible and can straightfowardly be adopted for the de- 
scription of other condensation phenomena such as for 
example ordinary Cooper pairing or a mixture of Cooper 
pair and a-particle condensation. 
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TABLE I. Comparison of the generator 


coordinate method calculations with 


experimental values. 


— nE a 


denotes the 



threshold energy for the decay into a-clusters, the values marked by * correspond to a refined mesh, see main text. 
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6.06 




o 3 + 


11.26 


2.6 


o 4 + 


12.05 


1.6xl0~ 3 


o s + 


14.0 


4.8 


o 6 + 


14.03 


2.0xl0 _1 



TABLE II. Observed excitation energies _B oxc and widths T of the lowest five + excited states in ie O 
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FIG. 1. Contour map of the energy surface E3 a (Ro,b) for 12 C. Numbers attached to the contour lines are the binding 
energies. 




FIG. 2. Contour map of the energy surface En a {Ro,b) for ie O. Numbers attached to the contour lines are the binding 
energies. 
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